From brittle to ductile fracture in disordered materials 
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We introduce a lattice model able to describe damage and yielding in heterogeneous materials ranging from 
brittle to ductile ones. Ductile fracture surfaces, obtained when the system breaks once the strain is completely 
localized, are shown to correspond to minimum energy surfaces. The similarity of the resulting fracture paths to 
the limits of brittle fracture or minimum energy surfaces is quantified. The model exhibits a smooth transition 
from brittleness to ductility, also present in how much plastic deformation is accumulated prior to fracture. The 
dynamics of yielding exhibits avalanches with a power-law distribution. 
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Understanding and characterizing the complexity in mate- 
rial failure is of great interest in both basic physics research 
and materials science applications. Brittle materials fail in the 
elastic regime, while ductile materials can locally accumulate 
plastic deformation prior to fracture and they are often able 
to withstand higher stresses before reaching fracture J2]. Nu- 
merical models of fracture often involve molecular dynam- 
ics simulations J2,[3t]. To go beyond these, brittle fracture 
with its simpler rheology has been extensively analyzed at 
the mesoscopic scale by both scalar and tensorial lattice mod- 
els J3]. This provides theoretical tools for understanding the 
intermittency observed in fracture- avalanches- and the ques- 
tion why crack surfaces roughen, to name some key prob- 
lems J3l . Other mesoscale models have been developed for 
plastic deformation and its accumulation in amorphous and 
disordered materials [5-9]. In plastic deformation one sees, in 
analogy to brittle fracture, general scaling laws like acoustic 
emission 1 13-12] and strain avalanche distributions 1 13- Tfih . 
A relevant phenomenon is the spatial localization of strain, 
sometimes noticeable as shear bands that can then induce fi- 
nal failure 1 17-2(J. One important issue in plasticity has been 
"optimization", as the yield surface may minimize the sum of 
certain local random quenched variables J5I E1T1 . 

In this Letter we study a scalar lattice model for fracture and 
plastic deformation in elasto-plastic heterogeneous materials. 
Depending on the material properties the model is able to ex- 
hibit fracture in situations that range from brittle to ductile be- 
havior. In analogy to dislocation dynamics and other scenar- 
ios of plastic deformation, for ductile systems the dynamics 
of yielding is characterized by strain bursts [ Io|-L5]. Finally, 
plastic deformation localizes into shear bands, until a crack 
develops. The model is a gene ralization of the well-known 
random fuse model (RFM) H22I1 for brittle fracture. Our duc- 
tile random fuse model (DRFM) is able to accumulate plastic 
deformation before complete failure. Depending on a model 
parameter, the model can interpolate between brittle failure or 
perfect plasticity depending on how ductile the system is, or 
how much it can yield. Due to the two intertwined dynam- 



ics of fracture and yielding the DRFM presents a very rich 
behavior, which we explore below. 

The model.- The RFM represents a scalar lattice electri- 
cal analog of the elasticity problem where the local stress 
(<7i), strain (e^), and local elastic modulus (Ei) are mapped 
to the current (Ij), potential drop (Vi), and local conductance 
(gi), respectively, of a network of fuses subject to an external 
voltage. Bus bar boundary conditions are imposed at the top 
and bottom of the system and periodic boundary conditions 
in the transverse directions. Each individual fuse i behaves 
ohmic (elastic), Ii = giVi (equivalently, Oi = Eiti), up to 
a local threshold current Ti, which is a uniformly distributed 
quenched random variable. 

The DRFM that we introduce here is based on the scalar 
tectonic model of Cowie et al. Q], but it could readily be ex- 
tended to tensorial models with more degrees of freedom such 
as beam models jfjj. Whenever a fuse reaches its threshold 
Ti, a permanent deformation is imposed to the element and it 
becomes elastic again. This defines a "healing cycle" of the 
individual element. In the electrical analogue this healing is 
done by imposing a voltage source (i.e. an electric battery) 
in series with the fuse to generate an opposite current through 
it so that elastic deformation is relaxed below threshold while 
plastic deformation accumulates in the element. The magni- 
tude of the local voltage source used to heal a fuse is linearly 
related to its threshold, Aj = f3Ti/gi where j3 is a param- 
eter that controls how much local deformation is allowed to 
be accumulated at each healing cycle. We repeat this heal- 
ing cycle a fixed number of times Nhc f° r eacn f use going 
above threshold (see top panel of Fig. Q]) until it definitely 
breaks. The number of healing cycles Nhc that each indi- 
vidual element can go through before failure is fixed for all 
the fuses and parametrizes in a simple and convenient manner 
the yielding characteristics (ductility) of the material. In our 
simulations we used /3 = 0.1 and the local elastic modulus 
or conductivity gi = 1, unless stated otherwise. Note that in 
the double limit Nhc 00 an d & — ► 0, one would obtain a 
elastic -perfectly plastic response. At each step of the Simula- 
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FIG. 1: (Color online) DRFM in a diamond lattice of size L = 
75 with disordered thresholds uniformly distributed in the interval 
[0.5, 1.5]. Top panel shows a typical stress-strain response curve for 
a single element of the lattice for j3 = 0.1. Left bottom panel shows 
global stress-strain curves for different typical disorder realizations 
with an increasing number of healing cycles and a decreasing median 
of the threshold distribution. For comparison, right bottom panel 
shows experimental data for amorphous steel alloys with increasing 
ductility from a to d (from Ref. 1241 ). 



tion both the fuse closest to threshold and the external voltage 
required to reach it can be exactly calculated. Therefore, the 
external voltage is increased exactly up to the point where the 
next fuse in the network reaches its threshold. The current re- 
distribution after the healing of a fuse can cause other fuses 
to also reach their thresholds. Therefore, avalanches of plas- 
tic events are observed similar to the strain bursts observed in 



experiments 11121 113L 12311 . 

Eventually fracture occurs as follows. At any given time in 
the simulation different sites of the network have gone through 
a different number of healing cycles, reflecting the spatially 
varying distribution of strain in the system. A fuse is forced to 
burn (break) after having gone through a fixed number of heal- 
ing cycles Nhc, so that this quantity parametrizes the capa- 
bility of the system to sustain local deformation. The fuse then 
irreversibly becomes an insulator and all the voltage sources 
that were imposed as plastic deformation are removed. When 
a continuous path of insulating bonds is formed, the system is 
disconnected and fails completely. 

Brittle versus ductile fracture surfaces.- We have studied 
the DRFM in two dimensions in a system of size L x x Ly. 
Bus bars at y = and y = L y + 1 impose an external voltage 
V across the system. Figure 1 (top panel) shows the stress- 
strain curve for an individual element of the network in the 
DRFM. Each stress drop corresponds to a healing cycle of the 
element and the amplitude of the zig-zag is proportional to j3. 
Global stress in this model is defined as a = I / L x and the 
strain is given by e = V/ L y , where / is the total current, V is 
the voltage drop, and L y = L x /2 + 1 = L + 1 for a system of 
size L in a "diamond" (45° tilted square) lattice. Figure[T|(left 
bottom panel) shows the global stress-strain curve for /3 = 
0.1 as the number of healing cycles Nhc before breakdown 



is increased and the median of the threshold distribution is 
progressively reduced. Also, in Fig. [T](right bottom panel) we 
plot stress-strain curves obtained by Gu et al. II24I1 in recent 
experiments with amorphous steel alloys of of Fe-Cr-Mo-P- 
C-B with different ductility produced by changing the metal- 
metalloid composition. The comparison illustrates how the 
index Nhc can parametrize in a simple manner the ductility 
of the experimental samples. Our results are also in excellent 
agreement with very recent experiments by Sun et al. iflrjll on 
ductile metallic glasses showing the cycles of sudden stress 
drop followed by elastic reloading associated with shear-band 
motions. 

Figure |2] (top panels) shows the resulting fracture paths in 
the DRFM for the same conditions and disorder configuration 
as in the simulations shown in Fig. Q] (bottom left panel). One 
immediately notices that the final fracture surface configura- 
tion depends on the number of healing cycles Nhc an d there- 
fore on the accumulated plastic strain. These surfaces are to be 
compared with the ones emerging from the perfect plasticity 
limit {Nhc - ► oo) 112 ill . As a reference, we plot the mini- 
mum energy (ME) surface l25ll and the perfectly plastic (PP) 
path found with the algorithm of Roux and Hansen H El for 
exactly the same disorder configuration. These two surfaces 
are known to minimize the sum of the local yield stresses and 
the stress flowing through the surface, respectively 12 ill . It 
becomes apparent that the deeper the system is allowed into 
the plastic steady-state (Nhc ^ 1)> the closer the resulting 
fracture surface is to the ME path for the same disorder config- 
uration. We define the overlap between two given paths {xi} 
and {yi} as (1/Z) J^- S(xi — yi), where Z is a normalization 
constant so that the overlap becomes unity for two identical 
surfaces. One can see in Fig. [2] (top panel) that the fracture 
surface tends to overlap with the ME surface as the ductility 
Nhc increases. Also, in Fig. [2] (bottom panel) we plot the 
overlap between the plastic fracture surface obtained with the 
directed ME surface for the same disorder realization. For 
the sake of comparison, we also compute the overlap with the 
corresponding purely brittle fracture surface fl22ll for the same 
disorder realization (i.e., setting Nhc = 0). The PP surface 
seems to be irrelevant for this problem and, consequently, the 
overlap is negligible for all Nhc ( n °t shown). 

Figure [2] demonstrates that the system is quasi-brittle for 
low values of Nhc where the fracture surface largely over- 
laps with the purely brittle fracture surface. In contrast, if the 
material is allowed to accumulate locally more strain (i.e. for 
larger Nhc) then the overlap with the brittle fracture rapidly 
decreases, while the path becomes progressively closer to the 
directed ME path. This is in agreement with the results re- 
ported in more complex models ]6[]. Moreover, the finite-L 
behavior of the overlap in Fig. [2] which decays as ~ log(l/L) 
with system size in both Nhc ~> oo (ME) and Nhc = 1 
(quasi-brittle) limits. 

Strain localization.- The dynamics of damage accumula- 
tion in the DRFM is summarized in Fig. [5] For quasi-brittle 
conditions, i.e. when we induce fracture after only a few heal- 
ing cycles (Nhc = 1 in the left top panel of Fig. |3J, the 




FIG. 2: (Color online) Top: ME, PP and DRFM fracture surfaces 
for the same disorder realization in a L = 75 system. One can see 
that the final fracture surface gets closer to ME as the system be- 
comes more ductile. Bottom: Average over 10 3 disorder realizations 
of the total spatial overlap of the DRFM final crack with both the 
corresponding ME and BF (Mhc = 0) surfaces for /3 = 0.1. 



damage in the form of broken bonds is spatially distributed 
throughout the sample in a randomly uniform fashion, in anal- 
ogy to the brittle RFM. In contrast, a ductile sample can ac- 
cumulate strain before any particular bond is broken. Right 
top panel in Fig. [3] shows how damage strongly localizes at 
the crack surface in a very ductile case with Mhc = 10 3 
healing cycles, which corresponds to an average accumulated 
plastic strain of order ~ MncPTiKgiLy) per site. In this 
case, the crack strongly localizes very close to the ME sur- 
face. Eventually, the damage localizes exactly at the ME path 
as Mhc — > oo (for a finite but large Mhc in a finite sample, 
cf. Fig. |2j). In the DRFM the yield localization behaves origi- 
nally randomly except that, due to stress enhancements, a de- 
gree of local "clustering" exists, and it increases slowly along 
the stress-strain curve until the final localization of damage 
and yielding as the fracture surface starts to be formed. Such 
a trend is analogous to what is seen in the brittle RFM for 
damage accumulation 0]. This is interesting since such local 
plastic strains could be measured in experiments. 

In the bottom panels of Fig. [3] we depict the spatial distri- 
butions of damage. The average damage profile p(y) is cal- 
culated from the fraction of broken bonds rib(y) along the 
y direction and is computed as p(y) = (ni,(y)) / L x , where 
the averaging is obtained by first shifting the damage profiles 
by the center of mass and then averaging over different sam- 
ples B26I1 . The damage becomes narrower for larger Mhc, as 
the material is allowed to accumulate more local irreversible 
strain. Prior to the growth of the final crack, there appears to 
be only local correlations in the damage similarly in the brit- 
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FIG. 3: (Color online) Top: DRFM fracture surface in a system of 
size L — 75 for Mhc = 1 (left) and Mhc = 10 3 (right) compared 
with the corresponding ME path for the same disorder realization. 
Bottom panels show the average damage profiles for the DRFM in a 
system of size L — 128 for different values of Mhc (left) and the 
corresponding data collapse for different system sizes in the case of 
large ductility Mhc = 5 x 10 3 (right). 



tie RFM M 

12711 . and the maxima in the profiles arise from 
the crack path. The bottom right panel of Fig. [3] shows the 
data collapse, p(y) = p(L/2)f[(y - L/2)/L Cme ], for the av- 
erage profiles of accumulated damage for very ductile sam- 
ples (Mhc = 5 x 10 3 ) with the local roughness exponent 
(me = 2/3 corresponding to ME. This quantifies how with 
an increased number of healing cycles fractures tend to only 
occur in the final crack surface, which progressively tends to 
the ME path. We also observe (not shown) that in the quasi- 
brittle limit damage is volume-like, as expected for RFM l26ll . 
while in the extremely ductile case damage scales as ~ L, as 
the failing elements are located at the one-dimensional crack. 

Avalanches of plastic events. - Indirect evidence by means 
of acoustic emission experiments | U 12|. and very recent di- 
rect micro-compression measurements Il3l4l5ll . have shown 
that strain bursts take place in the form of avalanches whose 
size distribution decays as a power-law with an exponent that 
lies within the range r = 1.4 — 1.6. On the other hand, the- 
oretical arguments 12811 and simulations in dislocation mod- 
els ificj I29I [30h and amorphous materials J3l suggest a uni- 
versal exponent r = 3/2. We measured strain avalanches 
as the total strain occurred in the system between two ex- 
ternal stress increments. In the DRFM the avalanche size 
is s = 5Zr=i A-i/Ly, which corresponds to the sum of the 
voltage sources Aj added between two external voltage incre- 
ments and n is the number of fuses involved in the avalanche. 
Figure[4](top panel) shows how the distribution of strain bursts 
for a sample of size L = 128 evolves towards the yielding 
point. As it can be observed the cutoff increases as the stress 
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FIG. 4: (Color online) Statistics of strain avalanches in the DRFM for 
j3 — 0.1 in a very ductile sample. Top: Avalanche distribution for a 
single realization in a system of size L — 128 collected up to differ- 
ent times in the evolution as marked in the inset. Bottom: Avalanche 
distribution in the plastic steady-state, well above the yielding point, 
for different system sizes averaged over realizations of disorder. In- 
set shows the corresponding data collapse. 



reaches the critical point, with an exponent that approaches 
t = 3/2. This is excellent agreement with recent measure- 
ments of the distribution of shear avalanches in ductile metal- 
lic glasses fjfill . On the other hand, well above the yielding 
point the material response is fully plastic and we measure a 
different strain avalanche distribution in this plastic steady- 
state, i.e. when the strain has already become completely 
localized and the average response in the global stress-strain 
curve in Fig.[T]is constant. To do this we start to record statis- 
tics of strain avalanches well above the yielding point and ob- 
tain V(s) = L- D (s/L D )- T f(s/L D ) with r = 1.1 ± 0.01 
and D = 0.2 ± 0.01, where / is a finite-size scaling function 
obtained by a data collapse for different system sizes (bottom 
panel, Fig. @). t —> 1 corresponds to the plastic flow regime 
where the path of yielding sites effectively separates the sys- 
tem into two parts. 

In summary, we have introduced a lattice model of elasto- 
plastic disordered materials which yield and finally fracture. 
The model exhibits a transition from purely brittle to fully 
ductile fracture, and includes both limits. A small accumu- 
lation of yield strain corresponds to quasi-brittle behavior. In 
contrast, as ductility is increased the resulting fracture paths 
gradually approach ME surfaces, and the damage decreases. 
The yielding process takes place in the form of avalanches of 
strain events that are found to be power-law distributed with 
an exponent r rs 1.5, which is to be compared with experi- 
ments reporting r = 1.4 — 1.6 jic 4l6ll and close to an earlier 
theoretical prediction, r = 3/2 [28}]. Interestingly, the DRFM 
also allows us to investigate avalanche dynamics in the plastic 



steady-state regime [a ~ O(e )] where we found r « 1.1, 
consistent with plastic flow and unbounded strain events in 
response to the external loading. Future directions from here 
would include studying the DRFM in three dimensions, and 
adding strain-hardening to the non-linearly elastic part of the 
element response. 
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